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INTRODUCTION:  One  of  the  basic  objectives  in  the  Interior  Ballistic 
Laboratory  of  BRL  is  the  development  of  an  overall  ballistic  model 
capable  of  simulating  the  performance  of  entire  weapon  systems  on 
the  computer.  Such  computer  simulations  would  be  capable  of  explor¬ 
ing  new  weapon  concepts  and  determining  the  pay  off  of  the  new  con¬ 
cept  prior  to  any  hardware  exploratory  development.  In  addition, 
guch  a  model  would  serve  to  guide  the  subsequent  experimental 
programs  needed  to  evaluate  the  feasibility  of  new  weapon  concepts. 

A  key  link  in  the  development  of  the  overall  ballistic  model  is 
the  solid  propellant  gun  gas  dynamics  computer  program.  This  program 
is  a  one-dimensional  Lagrangian  code  based  on  the  Richtmyer  Von 
Neuman  "q"  method;  but  unlike  the  basic  method  which  considered  the 
flow  of  gas  only,  the  present  program  considers  the  flow  of  a  fluid 
consisting  of  hot  gas  with  burning  propellant  grains  entrained  in 
the  gas. 

This  program  can  find  use  in  the  following  potential  areas  of 
gun  design: 

1.  Predicting  the  pressure  and  heat  transfer  distribution  in  a 
gun  where  the  propellant  chamber  diameter  is  greater  than  bore  dia¬ 
meter.  This  is  a  problem  to  the  gun  tube  designers  since  convention¬ 
al  interior  ballistic  theory  is  based  on  bore  diameter  chambers. 

With  the  emphasis  going  to  lighter  gun  tubes,  errors  in  pressure  and 
heat  transfer  distribution  results  cannot  be  compensated  for  by  gen¬ 
erous  safety  factors  in  gun  tube  thickness. 

2.  Predicting  the  acceleration  loadings  on  saboted  projectiles 
launched  from  high  velocity  guns.  High  velocity  launched  saboted 
kinetic  energy  penctrators  have  considerable  military  application  in 
antitank  warfare,  but  one  of  the  problems  in  their  use,  is  making 
the  sabot  light  enough  so  as  to  attain  the  desired  launch  velocity 
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and  yet  strong  enough  so  that  the  sabot  does  not  break  up  under  the 
high  acceleration  loadings  during  launch.  Prediction  of  the  acceler¬ 
ation  loading  on  these  projectiles  will  enable  the  sabot  designer  to 
design  light  weight  sabots  which  resist  the  stresses  encountered 
during  a  high  velocity  launch. 

BACKGROUND :  The  initial  study  on  the  gas  dynamics  problem  of  gun 
interior  ballistics  was  made  by  Lagrange  in  1793.  Lagrange  was 
responsible  for  formulating  the  Lagrange  problem  of  interior  bal¬ 
listics  which  referred  to  a  special  type  of  gun  in  which  the  projec¬ 
tile  is  prevented  from  moving  until  all  of  the  solid  propellant  in 
the  chamber  is  converted  to  a  gas.  The  problem  was  to  predict  the 
pressure  distribution  in  the  gas  after  the  projectile  began  to  move 
and  the  motion  of  the  projectile. 

Throughout  the  19th  century  and  the  early  part  of  the  20th 
century,  a  series  of  mathcmaticans  contribated  to  the  solution  of 
the  Lagrange  problem.  This  culminate  in  the  classic  work  of 
Love  and  Pidduck  (1).  In  this  work.  Love  c  irried  out  a  mathematical 
analysis  of  the  wave  motion  developed  in  the  propellant  gas  as  a 
result  of  the  motion  of  the  projectile.  Pidduck,  using  the  equatiens 
developed  by  Love,  computed  a  numerical  wave  solution  for  a  large 
caliber  gun.  Pidduck  and  later  Kent  (2)  found  a  limiting  solution 
to  the  Lagrange  problem  that  resulted  in  a  theory  that  the  ratio  of 
pressures  between  the  brcc-ch  face  and  projectile  base  was  held  to  be 
a  constant  during  the  firing  cycle.  This  ratio  is  a  function  of 
the  projectile  weight  to  propellant  weight  ratio  and  the  specific 
heat  ratio  of  the  propellant  gas.  This  theory  was  applied  by  Kent 
to  the  case  of  a  gun  in  which  the  propellant  continued  to  bum  while 
the  projectile  moved.  This  assumption,  however,  implied  that  the 
gas  density  in  the  propellant  chamber  at  start  of  projectile  motion 
was  non-uniform.  Since  the  experimental  results  of  interest  were 
maximum  pressure  and  muzzle  velocity;  this  assumption  was  sufficient 
for  most  gun  interior  ballistic  theories. 

With  the  advent  of  electronic  digital  and  analog  computers,  it 
became  possible  to  incorporate  these  interior  ballistic  theories 
into  computer  programs.  The  computer  programs  or  models  were  then 
capable  of  generating  gun  interior  ballistic  trajectories  which  then 
could  be  compared  with  measured  interior  ballistic  trajectories 
resulting  from  improved  gun  instrumentation.  Reference  (3)  describes 
such  a  computer  program  and  compares  predicted  interior  ballistic 
trajectories  with  corresponding  measured  values. 

While  the  computer  program  of  reference  (3)  is  able  to  simulate 
the  performance  of  any  type  of  solid  propellant  gun;  the  program  has 
certain  deficiencies.  One  of  these  deficiencies  is  that,  because  of 
application  of  the  Pidduck-Kent  constant  pressure  ratio,  the  pressure 
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on  the  projectile  base  during  the  projectile  engraving  is  too  low 
and  thus  predictions  of  the  effect  of  engraving  resistance  on  initial 
projectile  motion  in  the  tube  will  be  in  error.  Also  in  this  theory 
the  propellant  chamber  is  a  bore  diameter  extension  of  the  bore; 
thus  the  influence  of  chamber  dimensions  and  shape  on  gun  performance 
cannot  be  evaluated.  The  work  of  Seigel  (4),  while  confined  to  guns 
using  gas  propellant,  did  indicate  that  increasing  the  chambrage 
(propellant  chamber  diameter-bore  diameter  ratio)  would  increase 
the  muzzle  velocity  of  the  gun.  Another  deficiency  was  the  inabil¬ 
ity  of  the  program  to  predict  the  effect  of  change  in  gun  bore  cross- 
section  on  the  heat  transfer  distribution.  Experimental  work  at 
Cornell  (S)  had  indicated  that  the  heat  transfer  coefficient  varied 
from  breech  face  to  muzzle;  reaching  high  values  in  the  convergent 
cone  between  the  propellant  chamber  and  the  bore. 

For  the  above  reasons,  it  was  decided  to  develop  a  one-dimen¬ 
sional  gas  dynamic  gun  computer  model  which  would,  given  the  proper 
data,  be  able  to  predict  the  interior  ballistic  phenomena  discussed 
above . 

THEORETICAL  DEVELOPMENT:  The  theory  of  the  gas  dynamic  method  of 
gun  interior  ballistics  is  based  on  the  Lagrangian  form  of  the 
equations  of  time -dependent  fluid  flow.  In  general,  finite  differ¬ 
ence  methods  for  the  calculation  of  time- dependent  fluid  flows  have 
been  based  mostly  on  either  the  Eulerian  or  the  Lagrangian  form  of 
the  equations.  Although  these  forms  are  essentially  equivalent, 
the  Lagrangian  form  gives  more  information  (it  tells  where  each 
parcel  of  fluid  came  from  originally),  and  has  the  virtue  that  con¬ 
servation  of  mass  in  the  fluid  parcels  is  automatic  and  exact,  even 
in  the  finite  difference  approximations.  This  results  in  consider¬ 
ably  greater  accuracy  in  one-dimensional  problems.  In  addition, 
moving  boundaries  (such  as  the  base  of  a  projectile)  are  easy  to 
define  in  a  Lagrangian  system  and  are  more  difficult  to  define  in  an 
Eulerian  system. 

The  original  form  of  the  finite  difference  equations  describing 
one-dimensional  time-dependent  fluid  flow  were  derived  by  Richtmyer 
and  Von  Neuman  (6).  Since  in  a  solution  of  a  gas  dynamic  problem, 
shocks  are  likely  to  occur,  Richtmyer  and  Von  Neuman  devised  an 
artifical  dissipative  mechanism  of  such  form  and  strength,  that  the 
shock  transition  is  a  smooth  one,  extending  oveT  a  small  number  of 
space  intervals.  This  artifical  dissipative  mechanism  (artifical 
viscosity),  when  incorporated  in  the  finite  difference  equation, 
permits  all  dependent  variables  (pressure,  velocity,  etc.)  to  vary 
smoothly  through  the  shock  region  and  thus  suppress  any  numerical 
instabilities. 
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The  assumptions  used  in  the  model  are  as  follows: 

1.  Non-ideal  gas  obeying  the  Nobel-Abel  equation  of  state. 

2.  Frictional  losses  between  projectile  and  gun  bore  to  be 
either  a  constant  or  to  vary  with  projectile  travel. 

3.  Negligible  frictional  losses  between  propellant  gas  and  bore 
surface. 

4.  Heat  loss  to  be  a  function  of  gas  velocity  and  gas  temper¬ 
ature. 

5.  Unburned  solid  propellant  moves  as  fast  as  local  gas 
velocity. 

6.  Burning  rate  of  solid  propellant  is  a  function  of  the  local 
pressure  and  the  local  gas  velocity. 

The  following  Lagrangian  hydrodynamic  equations  are  used  in  the 
code: 

Momentum: 


3u 

at 


i(P-+-q) 

3M 


A(x) 


(1) 


Energy: 


These  equations  together  with  the  following  are  used  to  describe  the 
gas  dynamic  code: 

Nobel-Abel  Equation  of  State: 


P(V  -  b)  =  nRT 


(4) 


Propellant  Burning  Rate: 


B  P°  +  ky  v 


(5) 
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Propellant  Form  Function: 


Tu  develop  the  finite  difference  equations  used  in  the  computer 
model,  the  solid  propellant-gas  region  between  the  breech  end  of  pro¬ 
pellant  chamber  and  the  projectile  base  is  divided  into  J  segments  of 
equal  length.  This  is  illustrated  in  Figure  1.  Since  the  cross- 
sectional  area  A(x)  of  the  chamber  may  vary,  the  amount  of  solid  pro¬ 
pellant  in  each  segment  is  proportional  to  the  initial  volume  of  the 
segment.  Because  of  assumption  5,  the  total  mass  (solid  propellant, 
propellant  gas,  and  igniter  gas)  in  each  segment  will  remain  fixed 
throughout  the  remainder  of  the  problem. 


The  momentum  equation  in  finite  difference  form  is: 


and  the  gas  displacement  is: 


At  +  x 


(8) 


The  energy  equation  requires  a  special  procedure  to  place  it  in 
finite  difference  form,  because  of  the  presence  of  the  propellant 
chemical  energy  term  Q  and  the  desire  to  explicitly  solve  for  propel¬ 
lant  gas  temperature. 

In  finite  difference  form  the  energy  equation  is: 

A  (Q  -  E)  =  (P  +  q)  Av  +  AH^  (9) 


Using  the  procedure  in  reference  3,  we  define: 

QsJuz  C  dT  ♦  JmT  j  1  C  dT 
•o  v  1  'o  VI 


E  =  J 


m  zj  Cy  dT  +  Jmj  Cy 


(10) 
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Combining: 
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Q  -  E  =  J  m  z 


fTo  A 
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Q  -  E  =  J  m  z  Cy  (Tq  -  T)  +  jmj  Cy  (T^-T) 


Since : 
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,  and  JC 


V  -  (y-Dto  *  “Vj  (Yj-1)Tc 


then: 

O  p  _  miFZ  (T  -  T)  Vi  (T  -  T) 

Q  (Y-l)To  °  (Yj-l)^  °I 

4«  -  *>  ■  ,to  - 1*1)  -  (to  ■ r? 

*$nv#  "J-’TV 


r  ♦  pn  i 

(P  ♦  q)  4v  «  |(-j - J— ^  *  q"j  (V^1  -  v")  +  S«P  (I 

“l  '  C’  -  *  K}  0 

Combining  16,  17,  and  18  and  substituting  in  the  equation  of  state, 
we  arrive  at  the  following  equation  for  the  propellant  gas  tempera¬ 
ture: 

-n+1  £  -  *  -  »  n 

J  ~  I  1 


where: 
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(21) 
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H  =  (H"+1  -  H  n)  +  SH," 
3  Lj  J 
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(25) 


The  other  equations  are  converted  into  the  finite  difference  form  in 
a  similar  manner.  The  above  equations  are  discussed  in  greater 
detail  in  reference  8. 

A  computer  program  using  the  equations  of  the  theory  presented 
above,  was  written  for  the  two  ARDC  digital  computers  BRLESC1  and 
BRLESC2.  A  detailed  description  of  this  program  will  be  given  in  a 
subsequent  BRL  report  (9). 

INSTRUMENTED  CUN  SYSTEM:  In  order  to  verify  the  predictions  of  the 
gun  gas  dynamic  computer  program  described  above,  firings  were  under¬ 
taken  in  a  high  velocity  90-37mm  smooth  bore  gun.  The  dimensions  of 
the  gun,  the  projectile  weight,  and  characteristics  of  the  propellant 
are  given  in  Table  I.  Instrumentation  for  this  gun  is  illustrated  in 
Figure  2.  Pressure  instrumentation  on  the  gun  consists  of  pairs  of 
pressure  transducers  at  each  station  on  the  gun;  each  pair  consisting 
of  a  Minihat  strain  type  pressure  transducer  (10)  and  a  Kistler  Model 
607A  pressure  transducer.  Velocity  instrumentation  for  the  gun  con¬ 
sists  of  a  10  GHz  microwave  interferometer  to  measure  in-bore  projec¬ 
tile  motion  (11)  and  three  velocity  screens  in  front  of  the  gun  to 
obtain  muzzle  velocity.  Data  from  this  instrumentation  is  recorded 
on  analog  magnetic  tape.  This  data  is  subsequently  processed  on 
analog  to  digital  conversion  equipment  and  digital  computers,  using 
special  data  processing  routines,  to  produce  printed  digital  data. 

RESULTS  AND  DISCUSSION:  A  total  of  68  rounds  have  been  fired  in  the 
90-37mm  gun  to  date.  In  the  rounds  fired,  projectile  weight,  chamber 
volume,  and  propellant  weight  were  varied.  All  other  parameters  wei'e 
held  constant.  Muzzle  velocities  for  the  series  varied  form  4500  to 
7600  f/s  and  maximum  propellant  chamber  pressures  ranged  from  10,000 
psi  to  65,000  psi. 
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'At  the  present  time,  only  a  few  of  the  firings  in  this  program 
have  been  completely  reduced,  so  in  this  report,  only  the  results 
which  have  been  obtained  to  date  will  be  compared  with  predictions 
from  the  gas  dynamic  computer  program. 

To  obtain  predictions  from  the  gas  dynamic  computer  program  for 
the  rounds  which  will  be  discussed  here,  it  is  necessary  to  adjust 
certain  input  parameters  in  the  program,  values  for  which  are  not 
available,  so  that  predicted  maximum  breech  pressure  and  projectile 
muzzle  velocity  will  agree  with  the  corresponding  results  from  one 
of  the  experimental  rounds.  These  parameters  would  then  be  held 
constant  for  predictions  made  on  other  rounds  fired  in  the  gun.  For 
this  study  the  round  chosen  was  round  11.  The  constants  which  were 
adjusted  to  provide  the  agreement  between  theory  and  experiment  were 
engraving  resistance  profile,  heat  loss  coefficient  and  the  pro¬ 
pellant  erosion  constant  k  .  The  magnitude  of  the  resistance  pro¬ 
file,  and  the  other  two  constants  were  systematically  varied  until 
the  error  between  theory  and  experiment  for  maximum  breech  pressure 
and  muzzle  velocity  was  1%  or  less.  These  three  parameters  were 
then  held  constant  for  the  other  rounds  simulated.  The  values  of 
the  parameters  obtained  in  this  manner  are  given  in  Table  I.  No 
attempt  was  made  in  this  matching  process  to  obtain  agreement 
between  the  detailed  theoretical  and  experimental  interior  ballistic 
trajectories  other  than  maximum  breech  pressure  and  muzzle  velocity. 
Further  details  on  the  matching  process  are  given  in  reference  11. 

The  detailed  interior  ballistic  trajectory  of  the  round  used  in 
the  matching  study  is  illustrated  in  Figure  3.  In  the  figure, 
breech  face  pressure-time,  projectile  velocity-time,  acceleration 
pres sure -time,  and  projectile  displacement- time  are  plotted  both  for 
the  experimental  results  and  the  predicted  results.  The  experimen¬ 
tal  breech  pressure  curve  is  from  the  Kistler  pressure  transducer. 
The  velocity  and  projectile  travel  curves  are  from  the  interfer¬ 
ometer  data,  the  signal  of  which  faded  out  after  a  projectile  travel 
of  80  inches.  For  projectile  travel  beyond  the  80  inch  position, 
the  times  at  which  the  projectile  base  passed  under  the  pressure 
transducers  down  bore  were  used. 

To  correlate  the  theoretical  and  experimental  interior  ballis¬ 
tic  trajectories  in  time,  the  projectile  travel-time  curve  was  used 
as  a  base.  This  is  because  the  agreement  between  experimental  and 
theoretical  projectile  travel -time  curves  is  very  close  except  at 
start  of  motion,  which  is  experimentally  poorly  defined.  It  will  be 
noted  that  agreement  between  theory  and  experiment,  for  velocity- time 
is  good,  the  maximum  error  being  not  more  than  25%  between  the  two 
curves.  The  agreement  between  measured  and  theoretical  acceleration 
pressure- time  curves  is  also  good,  the  theoretical  curve  falling 
within  the  spread  of  the  experimental  data.  The  time  correlation 
between  the  measured  and  theoretical  breech  pressure-time  curves  is 
not  as  good.  At  maximum  pressure,  the  experimental  curve  lags  the 
theoretical  curve  by  about  0.5  ms. 
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The  reason  for  the  lag  between  the  two  curves  is  not  understood. 

A  probable  reason  may  be  that  the  gas  dynamic  computer  program  does 
not  correctly  simulate  gas  flow  and  propellant  burning  rate  in  the 
nozzle  connecting  the  propellant  chamber  and  the  gun  bore. 

Of  interest  to  gun  tube  designers  is  the  maximum  pressure  each 
portion  of  the  gun  tube  will  experience  during  the  interior  ballistic 
cycle.  Experimentally  this  is  determined  in  the  90-37mm  gun  by  deter¬ 
mining  the  maximum  pressure  recorded  from  each  of  the  pressure  trans¬ 
ducers  in  the  chamber  and  along  the  bore.  These  experimental  values 
are  plotted  against  pressure  transducer  position  referred  to  breech 
face  in  Figures  4, and  5.  These  plots  represent  rounds  9,  11,  and  12 
in  which  chamber  volume  projectile  weight,  and  propellant  weight  were 
varied.  It  will  be  noted  that  of  these  three  rounds,  one  used  a  44 
cubic  inch  chamber,  and  two  used  a  64  cubic  inch  chamber.  Propellant 
to  projectile  weight  ratios  varied  from  1.4  to  3.7  and  muzzle  veloc¬ 
ities  ranged  from  5000  to  S900  f/s.  Also  illustrated  in  these  plots 
arc  the  predicted  maximum  pressure  profile  from  the  gas  dynamic  code. 
In  general,  thz  theoretical  maximum  pressures  predicted  for  positions 
in  the  chamber  or  along  the  bore  either  fall  within  the  observed 
experimental  variation  in  maximum  pressures  or  fall  outside  of  the 
experimental  variation  by  an  error  of  r.ot  more  than  25%. 

In  table  II  are  given  projectile  weight,  propellant  weight, 
chamber  volume,  propellant  to  projectile  weight  ratio,  muzzle  veloc¬ 
ity,  and  maximum  acceleration.  Both  predicted  and  measured  muzzle 
velocity  and  maximum  acceleration  are  given.  The  percentage  error 
between  experimental  and  predicted  muzzle  velocities  is  not  greater 
than  5%,  and  the  error  between  experimental  and  predicted  maximum 
accelerations  is  not  greater  than  25%.  It  will  be  noted  that  pre¬ 
dicted  accelerations  are  higher  than  those  experimentally  observed. 

In  figure  6  theoretical  profiles  of  pressure,  gas  velocity,  heat 
flow  rate  to  bore  surface,  and  propellant  fraction  burnt  as  a  function 
of  position  in  the  gun  are  given  for  a  time  at  which  the  projectile 
has  moved  to  a  distance  of  33  inches  from  start  of  motion.  The  data 
is  taken  from  the  simulation  for  round  11.  Comparing  these  profiles 
to  the  gun  geometry  illustrated  at  the  bottom  of  the  figure,  it  can 
be  seen  that  the  30°  convergent  cone  between  the  propellant  chamber 
and  the  bore  creates  sharp  changes  in  the  slope  of  the  pressure  and 
velocity  profiles.  It  will  also  be  noted  that  the  heat  transfer  rate 
has  a  very  sharp  peak  at  a  position  just  down  stream  of  the  convergent 
cone.  This  is  the  area  in  the  gun  which  is  subjected  to  the  greatest 
amount  of  metal  erosion  and  according  to  reference  5  experiences  the 
greatest  amount  of  heat  addition.  The  fraction  of  propellant  burnt 
has  a  high  value  at  the  breech  end  and  a  lower  value  at  the  projectile 
base  end.  The  reason  for  this  is  that  propellant  in  the  breech  end 
is  subjected  to  higher  pressures  for  a  longer  time.  . 
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CONCLUSIONS: 

1.  The  gas  dynamic  computer  program  is  able  to  correctly  predict  the 
interior  ballistic  trajectory  (displacement,  velocity,  and  acceler¬ 
ation  versus  time)  of  the  projectile  while  in  the  gun,  provided  that 
data  for  resistance  profile,  erosion  coefficient,  and  heat  transfer 
coefficient  is  obtained  by  the  matching  of  an  experimental  round. 

The  velocity  agreement  is  with  an  error  of  10%  and  the  error  between 
measured  and  predicted  maximum  acceleration  is  less  than  25%. 

2.  Measured  breech  pressure-time  curves  show  a  time  displacement 
from  the  predicted  breech  pressure-time  curves.  This  displacement  is 
believed  to  be  due  to  the  possibility  that  the  one  dimensional  model 
does  not  correctly  simulate  the  two  dimensional  gas  dynamic  effects 
in  the  nozzle  between  the  propellant  chamber  and  the  gun  bore. 

3.  The  gas  dynamic  computer  program  is  able  to  predict  the  maximum 
pressures  to  which  the  gun  tube  is  subjected  along  its  length  within 
an  error  of  about  25%  of  maximum  pressure  at  a  position. 

4.  Heat  transfer  computations  using  the  gas  dynamic  computer  program 
indicate  that  the  point  of  maximum  heat  transfer  rate  is  very  close 
to  the  downstream  exit  of  the  nozzle  between  the  propellant  chamber 
and  the  bore. 

FUTURE  PLANS:  Plan.*  are  being  made  to  reprogram  the  gun  gas  dynamic 
program  from  the  FORAST  language  to  the  standard  FORTRAN  larguage  so 
that  the  program  can  be  used  by  other  agencies.  As  part  of  this 
reprogramming,  the  theory  will  be  revised  so  that  the  assumption  that 
unburned  propellant  moves  as  fast  as  the  gas  will  be  eliminated; 
motion  of  the  unbumed  propellant  will  be  governed  by  gas  drag  forces. 
Additional  changes  will  be  the  addition  of  provisions  for  the  burning 
of  more  than  one  propellant  in  the  gun  and  a  provision  for  simulating 
ignition  delays  of  portions  of  tve  propellant  bed. 

The  detailed  experimental  data  from  the  90-37mm  gun  and  other 
heavily  instrumented  guns  in  the  Interior  Ballistic  Laboratory  will 
be  used  to  verify  predictions  made  by  the  new  program  and  also  be 
used  to  modify  the  program  in  the  interest  of  obtaining  better  pre¬ 
dictions  of  new  and  unique  guns. 
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LIST  OF  SYMBOLS 
=  Viscous  resistance  constant 

*  Cross  sectional  area  of  j  gas  segment  interface  it 
position  x.. 

»  Propel lant-* gas  covolume 
«  Propellant  weight 

*  Igniter  Weight 

*  Specific  heat  at  constant  volume  of  propellant  gas 
a  Diameter  of  perforation  of  cylindrical  propellant  grain 
a  Outside  diameter  of  cylindrical  propellant  grain 
=  Internal  energy  of  propellant  gas 
=  Force  of  solid  propellant 
=  Force  of  igniter 
a  Gravitational  constant  386.09 
a  Heat  loss  to  bore  surface  from  j  gas  segment 
a  Mechanical  equivalent  of  heat 
=  Propellant  velocity  erosion  constant 
a  Length  of  propellant  grain 
=  Neight  of  propellant  gas  in  j  segment 
a  Neight  of  igniter  gas  in  j  segment 
=  Total  weight  of  j  segment 
a  Mass  propellant  gas 
a  Moles  of  gas 
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=  Number  of  perforations  in  cylindrical  propellant  grain 
=  Gas  pressure 
«  Pressure  in  j  gas  segment 

*  Chemical  energy  of  propellant 

=  Viscous  pressure  in  j  gas  segment 
=  Linear  burning  rate  of  solid  propellant 

*  Gas  constant  33572 

*  Time 

3  Initial  temperature  of  gas  ahead  of  projectile 

=  Gas  temperature 

«  Temperature  of  gas  in  j  segment 

»  Isoehloric  flame  temperature  of  propellant 

=  Gas  velocity 

=  Velocity  of  j  gas  segment 

=  Regression  distance  of  solid  propellant 

=  specific  volume  of  propellant  gas 

=  Volume  occupied  by  propellant  gas  in  j  segment 

*  Propellant  gas  volume 

■  Accumulated  work  done  on  or  by  j  segh.  n«. 

«  Displacement  of  j  gas  segment 

=  Weight  fraction  of  solid  propel Ian*  'urnt  ‘n  j  segment 
Propellant  burning  rate  exponent 
a  Propellant  burning  rate  coefficient 
=  Ratio  of  specific  heats  of  propellant  gas 
=  Integration  time  step  size 
=  Volume  of  segment 


TABLE  I 

Characteristics  of  90-37mm  Gun,  Projectiles,  and  Propellant 


Gun: 


Bore  Diameter: 

Maximum  PvojectiJe  Travel: 

Rated  Maximus  Pressure  of  Gun: 

Overall  Length  of  Gun: 

Chamber  Volumes  with  3  inch  I.D.  insert: 
Cone  Angle  between  Propellant  Chamber  and 
Projectile  Weights: 


1.456  in. 

239.4  in. 

100?000  psi 
22  ft-  8  in. 

44,  64,  and  84  cu.  in. 
Bore:  30° 

1/4  lb.,  1/2  lb.,  3/4  lb. 


Propellant: 

Type: 

Forte : 

^Specific  Heat  Ratio: 
Flame  Temperature: 
OoYolume: 

Web: 

Number  of  Fc^f orations: 
Igniter: 

Igniter  Weight: 


M2 

4,4ri3TP0  in- lb/lb 
1.222 
3372. BK 

29.  ©7  cu.  in. /lh 
.036  in. 

7 

Black  Powder 
.019  lb 


Matching  Parameters: 

Erosion  Constant  k  : 

Heat  Transfer  Coefficient  k. :  «> 

Resistance  Prefile:  *  ** 

Projectile  Travel,  in.  Resistance  Pressure,  psi 
0  800 

•4  4415 

•8  6989 

2*4  2483 

1.5  to  end  of  bore  700 


TABLE  II 

•  '•*  Firing  Data  for  90- 37mm 

Round  No. 

Projectile  Weight  -  lb. 

Propellant  Weight  -  lb. 

Propellant  -  Projectile  Weight  Ratio 
Chamber  Volume  in  3 
Muzzle  Velocity  Experimental  f/s 
Predicted  f/s 

Maximum  Acceleration  Experimental  kilo  g's 

Predicted  kilo  g*s 


Rounds 


Fig,  1,  Gas  fynamic  Model 


for  Solid  Propellant  Gun 
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Fife*  3»  Theoretical  and  Experimental  Interior  Ballistic 
Trajectory  for  Bound  11,  90-}7mm  Gun 


fig-  4.  Marimw  Pressure  Along  dun  Bore,  Bounds  9  and  12 


Tig,  6,  Ca3  Itynamic  Profiles  Along  90~37mm  Gun  Bore, 
Round  11  Theory  after  Projectile  Travel  of 
33  inches. 


